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Abstract 

Ricci and sectional curvatures of twisted flux tubes in Riemannian manifold are computed 
to investigate the stability of the tubes. The geodesic equations are used to show that in 
the case of thick tubes, the curvature of planar (Frenet torsion-free) tubes have the effect 
ct of damping the flow speed along the tube. Stability of geodesic flows in the Riemannian 
twisted thin tubes (almost filaments), against constant radial perturbations is investigated by 
using the method of negative sectional curvature for unstable flows. No special form of the 
flow like Beltrami flows is admitted, and the proof is general for the case of thin tubes. It is 
found that for positive perturbations and angular speed of the flow, instability is achieved , 
since the sectional Ricci curvature of the twisted tube metric is negative. PACS numbers: 

02.40.Hw:Riemannian geometries 
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I Introduction 



The stability of geodesic flows have been recently investigated by Kambe [1] by making use of 
the technique of Ricci sectional curvature [2], where the negative sectional curvature indicates 
instability of the flow, while positivity or nul indicates stability. In the case of instability the 
geodesies deviate from the perturbation of the fluid. Following the work of D. Anosov [3] on 
the perturbation in geodesic flows in three-dimensional Riemannian geometry, in this paper 
the sectional Riemann curvature of the geodesic flow for a Riemannian flux tube [4, 5], where 
the axis of the tube flow possesses Frenet curvature and torsion. In the approximation of a 
thin tube where the radius of the tube is almost null, we show that the flows are unstable, 
against orthogonal perturbations, which is equivalently due to the negativity of the sectional 
Ricci sectional curvature. Throughout the paper the ellegant coordinate-free language of 
differential geometry [6] is used. The paper is organized as follows: Section II presents a 
brief review of Riemannian geometry in the coordinate free language. Section III presents the 
geodesic flow computation of the Christoffel symbols for the thick flux tube, where we show 
that the curvature of the tube axihe speed of the flow. Section IV presents the computation 
of the instability of Riemannian tube flow. Section V presents the conclusions. 

II Ricci and sectional Riemann curvatures 

In this section we make a brief review of the differential geometry of surfaces in coordinate-free 
language. The Riemann curvature is defined by 



where XeTM. is the vector representation which is defined on the tangent space TAi to the 
manifold M. Here VxY represents the covariant derivative given by 



which for the physicists is intuitive, since we are saying that we are performing derivative along 
the X direction. The expression [X, Y] represents the commutator, which on a vector basis 



R(X, Y)Z := V X V Y Z - V Y V X Z - V [X , Y] Z 



(11.1) 



V X Y = (X.V)Y 



(11.2) 
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frame e.\ in this tangent sub-manifold defined by 

X = X k e k (II.3) 

or in the dual basis d k 

X = X k d k (II.4) 

can be expressed as 

[X,Y] = (X,Y) k d k (II.5) 
In this same coordinate basis now we are able to write the curvature expression (1 1.1) as 

R(X, Y)Z := [R l jkp & X^di (II.6) 

where the Einstein summation convention of tensor calculus is used. The expression R(X, Y)Y 
which we shall compute bellow is called Ricci curvature. The sectional curvature which is very 
useful in future computations is defined by 

where S(X,Y) is defined by 

S(X,Y) := ||X|| 2 ||r|| 2 - < X,Y > 2 (11.8) 
where the symbol <, > implies internal product. 
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Ill Geodesic equations in Riemannian tube metric 

In this section we shall consider the twisted flux tube Riemann metric. The metric g(X,Y) 
line element can be defined as [4, 5] 

ds 2 = dr 2 + r 2 d6 R 2 + K 2 {s)ds 2 (III.9) 

This line element was used previously by Ricca [4] and the author [5] as a magnetic flux tubes 
with applications in solar and plasma astrophysics. This is a Riemannian line element 

ds 2 = g ij dx i dx j (III. 10) 

if the tube coordinates are (r, 6 R , s) [4] where 9(s) = 9 R — J rds where r is the Frenet torsion 
of the tube axis and K(s) is given by 

K 2 (s) = [1 - rK(s)cos6(s)} 2 (III.ll) 

Let us now compute the geodesic equations 

^ + ry/ = o (in.12) 

where v s = ^ and v e = ^ the Riemann-Christoffel symbols are given by 

r l ifc = \g il [9i3,k + 9ikj - 9jk,i] (in. 13) 

The only nonvanishing components of the Christoffel symbols or Levi-Civita [6] connections 
of the flux tube are 

r 2 21 = - (IH.14) 

r 

r 2 33 = - K ^ Ksme (in. is) 

r 

r 3 n = ^-[rsinO + -cosO] (III. 16) 

2K k 

r 3 33 = K- X K' = -r 3 n (111.17) 

A simple example can be given by writing the geodesic equation for the untwisted tube, where 

v e = 0, as 

s + r 3 n [r 2 - s 2 } = (111.18) 



where substitution of the component r 3 n above yields 

dlnv s = -d(lnn) (III. 19) 

where we use the chain rule of differential calculus, = — . Solution of equation (111.19) 

Vs Oil Vs 

is 

Vs = v^k- 1 (111.20) 

where v is an integration constant. This solution tells us that when the curvature tends to 
oo the velocity along the tube axis vanishes. Physically this means that the curvature acts as 
a damping to the flow, along the tube. In the next section we investigate in some detail the 
stability of an incompressible or volume preserving flow, using the method of the sign of the 
Ricci sectional curvature. 



IV Ricci sectional curvature and flow stability 

One of the most important features of the investigation of the stability of flows in the Euclidean 
manifold £ 3 , is the comprehension of the fact that the covariant derivative in the flow curved 
manifold is given by the gradient operator in curvilinear coordinates. Thus to compute the 
Ricci sectional curvature above, we need to make use of the grad operator in the twisted flux 
tube Riemannian metric given by 

V=[d r ,r- 1 d dR ,K- 1 d s ] (IV.21) 

Since the axis of the tube undergoes torsion and curvature, we need some dynamical relations 
from vector analysis and differential geometry of curves [7] such as the Frenet frame (t,n,b) 
equations 

? = kvl (IV.22) 

n' = -Kt + rb (IV.23) 

b' = -rn (IV.24) 

and the other frame vectors are 

e r = ncosO + bsinO (IV.25) 
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eg = —nsinO + bcosO (IV.26) 
d d e e = -n[(l + T- X K)sinQ + cosO] - b\cos6 + smO] (IV.27) 
Let the constant perturbation be given by 

X = u r l e r (IV.28) 

The upper index one in this expression refers to the fact that the background original value of 
u r was considered as u r ° = to form the tube. The other variable Y is given by 

Y = u e e + u s t (IV.29) 

Therefore to compute the Ricci tensor step by step we start by the term 

V X Y = u r (1) d r [u e ee + u s ^ (IV.30) 

which vanishes since we addopt here the approximation u r ^d r [ue] ~ along with the same 
relation to the radial partial derivative of u r . So 

VyVxY w (IV.31) 

Now the second term in the Ricci tensor is 

V x VyF = -u r (1) r- 2 u e [ (yl + T ^ (ncosO + bsin9)\ (IV.32) 

T 

where we have used the approximation of the thin tube where K(s) ~ 1 and r w 

[X.Y] = u r {l) [r- l u e - tu s ] [e r - (IV.33) 

which implies that 

V[x,Y]Y = « r (1) [r^ug - TU S ] [-r^ugd^g + Kn] (IV.34) 

The Ricci tensor is 

R(X, Y)Y = -u r {1] r^uer' 1 d s e e (IV.35) 
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In the previous computations we have made use of the imcompressibility of the flow 

V.m = (IV.36) 

which is 

d s ue = ttkuq ps (IV. 37) 

since r ps on the RHS of equation (IV. 37). The sectional curvature is thus 

<R(X,Y)Y,X: u e [l + r(s)Kcos9] 
K{X > Y) ~ S{XX) " " ruW r W + u s V 



when the tube is strongly twisted, uq » u s 2 thus the sectional Ricci curvature is 

< R(X,Y)Y,X> [1 + t(s)kcos9] 
Y) = S(X^) = ruV r u e (IV - 39) 

when the tube, besides is planar or torsion vanishes the last expression reduces to 

< R(X,Y)Y,X> 1 rjV4n , 

K(X ' Y) = S(X-Y) = -^Ue (IV - 40) 

Note that when both angular velocity and perturbation both keep the same sign, the sectional 
curvature K(X, Y) is negative and the flow along the Riemannian flux tube is unstable. There 
is singularity in this sectional Riemannian curvature in r ps 0. 



V Conclusions 

An important issue in plasma astrophysics as well as in fluid mechanics is to know when a fluid, 
charged or not, is unstable or not. In this paper we discuss and present incompressible flows 
and investigate their stability. Instability is obtained even before the singularity is achieved. 
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